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MINIMUM VARIANCE UNBIASED AND MOMENT ESTIMATORS OF 



CARRIER FREQUENCY OFFSET IN MULTI-CARRIER SYSTEMS 



CROSS-REFERENCE TO RELATED APPLICATION 



5 



This application is related to a co-pending application entitled "Globally 



Optimum Maximum Likelihood Estimation of Joint Carrier Frequency Offset and 
Symbol Timing Error," U.S. Ser. No. 09/496,890, filed on March 2, 2000 5 assigned to the 
assignee of the instant invention, and the disclosure therein is hereby incorporated by 
reference into the instant application. 



1 . Field of the Invention 

The present invention relates to a method for data communication of signals in 
units of a frame using an orthogonal frequency division modulation (OFDM) algorithm, 
15 and in particular, to methods of estimating carrier frequency offset error at the receiver. 

2. Description of the Prior Art 

OFDM system is a viable modulation scheme for data transmission time varying 
dynamic channels. However, it is known that performance of such system is highly 
susceptible to non-ideal synchronization parameters. Specifically, symbol timing and 
20 carrier frequency offset become an increasingly important issue in implementation of 
OFDM systems for practical applications. It is known that carrier frequency offset 
deteriorates performance of OFDM systems by introducing interference among the sub- 
channels. To overcome this imperfection, various compensation methods for estimation 
and correction of synchronization parameters have been proposed. In order to compare 
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the performance of these estimators, it is required to define a single number representing 
the goodness of the estimate. Assuming that all estimators are unbiased, i.e., expectation 
of the estimate is equal to the parameter, the variance of the estimator is used as a global 
measure for performance comparison of these estimators. 
5 Cramer-Rao lower bound (CRLB) is a fundamental lower bound on the variance 

of the estimators and the unbiased estimator whose variance equals CRLB is called 
efficient. When the evaluation of efficient estimator is not possible, it is desirable to 
obtain an estimator in which its performance becomes as close as possible to the CRLB 
fundamental bound. The estimator which is closest in performance to the CRLB 

10 estimator is known as a minimum variance unbiased (MVU) estimator. 

Categorically, the previously proposed methods for synchronization of OFDM 
systems can be classified into two main subclasses, namely minimum mean square error 
(MMSE) and maximum likelihood (ML) estimators. In MMSE approach, the estimator 
uses the information provided by the reference signal (pilot tones) in order to minimize a 

15 cost function associated with the synchronization parameters. A salient feature of this 
approach is that no probabilistic assumptions are made with regard to the data. Although 
MMSE estimators usually result in a tractable (globally stable) and easy to implement 
realization, no optimal criteria (probabilistic) is associated with these estimators. Also, 
since part of the transmitted information is allocated to the reference pilots, the 

20 bandwidth efficiency of these methods is lower in comparison to the non-pilot schemes. 

On the other hand, ML estimators provide the estimate of the unknown parameter 
subject to minimum probability of error criteria. Although not exactly efficient, ML 
estimators are asymptotically MVU, i.e., their variance attains that of MVU estimator as 
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the length of data record goes to infinity. However, due to the physical constraints, 
systems with infinitely long data records are not feasible for implementation purposes. 

P.fL Moose, in "A Technique for Orthogonal Frequency Division Multiplexing 
Frequency Offset Correction," in IEEE Trans. On Communications, Vol 42, No. 10, pp. 
5 2908-291 3, Oct. 1 994, describes the use of a retransmission technique in order to reveal 
the frequency offset parameter in the likelihood function of the received signal. Due to 
the redundancy introduced by repeating the data block, the data rate efficiency is 
decreased by a factor of two. To avoid this imperfection, a ML estimator based on cyclic 
prefix (CP) is described by J. van de Beck, M. Sandel and P. O. Borjesson, in "ML 

10 Estimation of Timing and Frequency Offset in OFDM Systems," IEEE Trans. On Signal 
Processing, Vol. 45, No. 3, pp. 1800-1805, July 1997. In this approach, the side 
information provided by the CP is used to obtain the likelihood function for joint 
estimation of symbol timing error and frequency offset in an OFDM system. 

The likelihood function described in the Moose reference does not globally 

15 characterize the observation vector over the entire range of the timing offset. 

Consequently, the ML estimator proposed based on this likelihood function would result 
in considerable performance loss over a finite range of timing offset interval. 

Currently, there is increasing interest in multi-carrier modulation (MCM) for 
dividing a communication channel into several subchannels and transmitting many 

20 subcarriers through a single band using frequency division multiplexing (FDM) 
techniques. In the MCM method, however, because several subcarriers occupying a 
narrow frequency domain are transmitted at one time, a relatively longer symbol period 
results compared with a single carrier modulation method. The MCM method has, owing 
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to such characteristics, the advantages that equalization is easily performed and that it has 
immunity to impulse noise. OFDM is a type of the MCM designed to maximize the 
working frequency efficiency by securing orthogonality among the multiplexed 
subcarriers. OFDM is applied to mobile radio channels to attenuate multipath fading. 



parallel data are carried out using the Fast Fourier Transform (FFT). It is required that the 
sampled data be sent in predetermined frames, having passed through a FFT routine, been 
time-division multiplexed, and transmitted, then restored at the receiving end. However, 
if the synchronization is in error in the course of restoring the frame, the signals 
10 demodulated after the FFT will be influenced by interchannel and intersymbol 
interference. Accordingly, the problem of synchronization in reforming the frame, 
especially any joint carrier frequency offset or symbol timing error, must be addressed as 
a matter of importance. 



15 that the process of synchronization is not only very complex, but the synchronization is 
not realized rapidly. 



Motivated by the sub-optimum performance of estimators of the prior art, a 
likelihood function for joint estimation of carrier frequency offset and symbol timing 
20 error of OFDM systems is disclosed in the related, copending application; and a new 
optimum ML joint estimator is disclosed therein. In order to reduce the variance of that 
ML estimator, a new class of MVU estimators for frequency offset estimation of PFDM 
systems is disclosed herein. 



5 



In an OFDM transmitting/receiving system, modulation and demodulation of 



Conventional synchronization methods as above-described encounter problems in 
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There is disclosed to exist but one function of sufficient statistic which provides 
the MVU estimate of the frequency offset. The estimator provided by the instant 
invention is a closed form expression; providing an estimator which is a function of data 
statistic. Consequently, it does not suffer from converging to multiple local minima, a 
5 problem which arises in ML technique with nonconvex loglikelihood functions. 

The advantages of the instant MVU estimator over the class of previously 
proposed estimators are two; first, it is MVU, therefore its variance is minimum among 
the entire class of estimators which use the same probabilistic measure. Secondly, it 
provides a closed form expression for mapping the statistics into the estimation domain. 
10 The former property assures optimality of the estimator, while the later facilitate the 
closed loop analysis of the system. 

Accordingly, the present invention is directed at a synchronization method that 
substantially obviates on or more of the problems due to limitations and disadvantages of 
the prior art. To achieve these and other advantages, and in accordance with the purpose 
15 of the invention as embodied and broadly described, there is provided a method of 
estimating carrier frequency offset error in a received sample bit stream including an 
observation vector (OV), having an observed carrier frequency timing offset e, and a 
plurality of data-symbol frames, having a symbol timing offset error d. The method 
comprises the steps of generating a probability density function (PDF) based on the OV 
20 and generating from the PDF an estimate of carrier frequency offset error, 8mvu I & being 
a minimum variance unbiased (MVU) estimator. 

In another aspect of the invention, the OV comprises an L-bit cyclic extension 
portion and a first and a second N-bit synchronization frame, and wherein the PDF 
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10 



15 



comprises a first term, pi, based on the timing offset & being within the span 1 to N and a 
second term, p2, based on said timing offset 3 being within the span N+l to N+L. In yet 
another aspect of the invention the received bit stream has uncorrelated independent 
identically distributed random signal and noise sequence variables with power of <y] and 

cr 2 n , respectively, wherein the OV is denoted x, and wherein the MVU estimator smvu I s 
is the conditional expectation of a second moment estimator, said second moment 
estimator given by 



1 & + L ~ l r i 
1 2 x[k]x*[k+N] 



Lcr} k = 3 



In another aspect of the invention the MVU estimator Smvu I a is given by 





r 




MVU|,9 = E 


e 





Yx[kle*[k + N] 
La 1 ^ 



= — ln£ 



's k=9 



7j(x,5) | 



In Zcr 2 



where 3 is the imaginary operator and where 

S-N-l MIA 

Y,ik}cik+N\+ £#*[*+^ N+l<$<N+L 



i;(x,5)= 



The invention also provides a method of synchronizing a received sample bit 
stream, comprising the steps of transmitting at a transmitter the bit stream including an 
observation vector (OV), receiving and sampling at a receiver the bit stream, the sampled 
bit stream including, the OV with an observed carrier frequency offset s, and a plurality 
of data-symbol frames, having a symbol timing offset error S; generating a probability 
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density function (PDF) based on the OV; generating from the PDF an estimate of carrier 
frequency offset error, smvu |s> being a minimum variance unbiased (MVU) estimator; 
and synchronizing the received bit stream by the MVU estimate of carrier frequency 
offset. 

5 Further aspects of the synchronization method of the invention include wherein 

the OV comprises an L-bit cyclic extension portion and a first and a second N-bit 
synchronization frame, and wherein the PDF comprises a first term, pi, based on the 
observed timing offset d being within the span 1 to N and a second term, p2, based on the 
observed timing offset S being within the span N+l to N+L. In yet another aspect of the 

10 inventuion, the synchronization method provides that the received bit stream has 
uncorrelated independent identically distributed random signal and noise sequence 

variables with power of a] and a\ , respectively, wherein the OV is denoted jc, and 
wherein the MVU estimator smvu|s is the conditional expectation of a second moment 
estimator, the second moment estimator given by 



15 e=^-\?A 

271 



j 3 + L-l 



Lcrj k = 3 



£ x[k]x*[k + N] 



Further aspect of the invention, the synchronization method wherein said MVU 
estimator Smvu | a is given by 



£ mvu 1 ,9 = E 







£ 

V 





271 \La s k=s 



f , 9+L-l 



In Lai 



20 where 3 is the imaginary operator and where 
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k=0 k=S 



In another aspect of the invention, A method of estimating carrier frequency 
offset error is provided for a received sample bit stream including an observation vector 
(OV) ? having an observed carrier frequency timing offset s, and a plurality of data- 
5 symbol frames, the OV comprising an L-bit cyclic extension portion and a first and a 
second N-bit synchronization frame, comprising the steps of generating the expected 
value of the autocorrelation of the Mi entry of the OV; and generating from the expected 
value an estimate of carrier frequency offset error, e mom , being a moment estimator. 

In yet another aspect of the invention, the moment carrier frequency offset error 
10 estimation method further includes wherein the OV is denoted x, and wherein the 

moment estimator s mom is given by s mom - {inT^ (x) 

where the statistic (x) is defined as 

t,(x) = ^^ytvr*+jvi- 

3 La} yt=0 L J L J 

A method of synchronizing a received sample bit stream, is also provided by the 
15 invention comprising the steps of: transmitting at a transmitter said bit stream including 
an observation vector (OV); receiving and sampling at a receiver said bit stream, said 
sampled bit stream including, said OV with an observed carrier frequency offset s, and a 
plurality of data-symbol frames, having a symbol timing offset error 8 ; 
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generating the expected value of the autocorrelation of the kth entry of the OV; 
generating from the expected value an estimate of carrier frequency offset error, s m om? 
being a moment estimator; and synchronizing the received bit stream by the moment 
estimate of carrier frequency offset. 

Further aspect of the invention provides a synchronization method wherein the 
OV is denoted x, and wherein the moment estimator e mom is given by 



In yet another aspect of the invention, a method of estimating carrier frequency 
offset error in a received sample bit stream including an observation vector (OV), having 
an observed carrier frequency timing offset s ? and a plurality of data-symbol frames, the 
0 V comprising an L-bit cyclic extension portion and a first and a second N-bit 
synchronization frame, comprising the steps of: selecting one estimator method from a 
group of three: a maximum likelihood (ML) estimator, a minimum variance unbiased 
(MVU) estimator, a moment estimator, based on said OV; and generating from the 
selected estimator an estimate of carrier frequency offset error. 

Additional features and advantages of the invention will be set forth in the 
description which follows, and in part will be apparent from the description, or may be 
learned by practice of the invention. The objectives and other advantages of the invention 
will be realized and attained by the methods particularly pointed out in the written 
description and claims hereof, as well as the appended drawings. 




where the statistic T^(x) is defined as 
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BRIEF DESCRIPTION OF THE DRAWINGS 
FIG, 1 is a signal flow block diagram depicting the frequency recovery loop for 
5 estimation of frequency offset. 

FIG. 2 is a block diagram depiction of the ML estimator disclosed in the cross- 
referenced, copending application and the MVU and moment estimators of the instant 
invention. 

FIG. 3 is a graph comparing the computer simulations of the prior art ML 

10 frequency offset estimation and the ML frequency offset estimation disclosed in the 

related, copending application. 

FIG. 4 is a graph comparing the computer simulations of the MVU frequency 

offset estimator of the instant invention with that of the CRLB estimator. 

FIG. 5 graphically depicts the closed loop performance of the MVU estimator of 

15 the instant invention by comparing the frequency offset estimate of the closed loop MVU 

estimator against an analytical derivation thereof 

DETAILED DESCRIPTION OF THE INVENTION 

The cross-referenced co-pending application entitled "Globally Optimum 

Maximum Likelihood Estimation of Joint Carrier Frequency Offset and Symbol Timing 

20 Error," U.S. Ser. No. 09/496,890, filed on March 2, 2000, incorporated herein by 

reference, discloses a probability density function (PDF),p, which globally characterizes 

an observation vector x according to the equation 

p(x 9 e, ff) = a (x, s 9 9){U[S - 1] - U\& - N - 1]) 

+ p 2 (x,s,&)(U[S-N-l]-U[S-N-L + l]) K ) 
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where e and 9 are the frequency offset and symbol timing error introduced by the 
synchronization mismatch in the carrier frequency and symbol timing, respectively. And 
Pi and p 2 are probability density functions derived and described in the related copending 
application. 

Further derived and described therein are the statistics Ti and T 2 



T, (x^- 



and 



^4fc]x*[A:+jV] 1<9<N 

MA 



(15) 



2 \x[k]\ 2 + \x[k + N]\ 2 1< & < N 
2>[*]| + + Z \x[k]\ 2 +\x[k + N]\ 2 N + l<S<N + L 

NandL being the number of samples per frame and the number of cyclically-extended 
samples appended at the beginning of each frame, respectively, all as further described in 
the related, copending application. 

The MVU estimator of the instant invention is derived by resorting to the theory 
of sufficient statistics. The first step in deriving the MVU estimator is to obtain the 
sufficient statistic for the PDF given in (14). The sufficient statistic is known to be a 
function of the observation vector, namely 7(x), such that the conditional PDF of the 
observation vector given T(x) does not depend on the unknown estimation parameters 
[£,&]. Evaluating the sufficient statistic is a formidable task for the broad class of PDFs, 
however the Neyman-Fisher Factorization theorem can be used for identifying the 
potential sufficient statistics. According to this theorem, if the PDF can be factored in the 

n 
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form g(T ( x), s , 9, )/?(x) where g is a function depending on x only through T(x) and h(x) 
is a function depending only on x, then T(x) is a sufficient statistic for estimation of the 
parameters s and 9 . By reformulating the PDF given in 14 to 

25Rp 2 ^r 1 (x,5y-H 2 r 2 (x,5) 



5 there is a direct dependency between the parameter 9 and the statistics Ti( x, 9 ) and 
T 2 (x, 9) . Based on this observation, the Neyman-Fisher theorem fails to provide a 
sufficient statistic for estimation of 9 , we can factor the PDF into 

5Rp' 27K r i (x,5)|a| 



H«l 2 f °-h°-%) 

p(x,e\ s )=e ^ A > A 2 (x) (17) 
Clearly then, 7;(x,«9) forms a sufficient statistic for estimation of the parameter e . 
1 0 Next, application of the Rao-Blackwell-Lehmann-Scheffe (RBLS) theorem is 

used to find the MVU estimator. According to this theorem, if s is an unbiased 
estimator of s and T(x) is a sufficient complete statistic for s then s = E(s\ m ) is a 
valid, unbiased, MVU estimator of s . 

In applying the above theorem, we need to obtain an unbiased estimator of 
1 5 s , termed s , and determine the conditional expectation of this estimator given the 
statistic T x (x, An appropriate candidate for the unbiased estimator of s can be 

obtained from the statistical moments of the random vector x. The autocorrelation 
function of the observation vector x, E, yields that the second moment of the random 
variable x[k] with keCl, satisfies the following identity 
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E{x[k\c*[k+N]}=ah~ j27l£ (18) 
Having this observation, one uses the second moment estimator as an unbiased estimator 
for s as given by 



2n 



Ley 



£ x[k]x*[k+N] 



s k = & 



(19) 



5 where 3 is the imaginary operator. In deriving the above estimator, £(x[A:Jc * [k + iV"]) 
was replaced by its natural estimator jYZT^ *l k+N \ !t is straightforward to verify 
that this estimator is unbiased as it satisfies the condition 



E{e} = —\nE- 
In 



S+L-l 



La 2 



(20) 



Next, one obtains the conditional expectation of s given the sufficient statistic r i(x s) as 
10 follows 





f 




MVUj^ = E 


£ 







's lc=S 



i 3 W 

2n La 2 . 



(21) 



It is important to emphasize that since the underlying PDF given in (14) belongs 
to the exponential family of PDFs, then the sufficient statistics T x {^,&) forms a complete 
1 5 statistic for estimation of the parameter e . Therefore, the mapping function obtained 
from applying RBLS theorem, namely lnr(x,5), is but one function of the statistic 
^(x, &) and no other estimator with the same statistic can result in a lower variance with 
respect to MVU estimator. 
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A, Cramer-Rao Lower Bound 

It is known that under broad conditions, the variance of any unbiased estimator of 
a nonrandom parameter s satisfies the CRLB as 



™^MVU\&)>^ 



to 



(22) 



5 where l(e) is the Fisher Information given by 



1(e) = -E 



d 2 \apr(x,£\ 3 ) 



Be 2 



(23) 



Substituting (17) into (23), after some algebraic manipulations, the CRLB of the MVU 
estimator becomes 



var 



fc )_ (i-i«rk^ 2 ) 

{mvuiJ (2^) 2 |a| J e:{r 1 (x„9)}" 



i+ 



l 



SNR 



-1 



(2x) 2 L 



(24) 



A (J o 

10 where SNR = — ^- is the signal to noise ratio at the receiving end. 

B. Closed Loop Performance 

FIG. 1 is a signal flow block diagram depicting the frequency recovery loop for 
estimation of frequency offset according to the instant invention. A closed loop system is 
obtained by feeding back the information obtained from the estimator into a sampler 
15 block (bootstrap). A sampler 10 updates its frequency at the beginning of each 

observation vector (every (N+L) samples). The sampler 10 receives the latest signal and 
produces the sampled symbol stream x[k] therefrom which is conducted to a block 12 and 
a block 14 performing the autocorrelation operation on the observation vector as required 
by equation (19). The result of the autocorrelation operation is conducted to a moving 
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average (MA) filter block 16 and therefrom to a block 18 performing the natural 
logarithm function and thence to a block 20 performing the imaginary operation function, 
as required by equation (19). To match the various sampling frequencies used in the 
system, a down sampler (decimator) block 22 is used prior to return to the sampler 10 to 
5 produce the minimum variance unbiased estimator, s [m]. Finally, a gain block (G) 24 is 
used to control the closed loop characteristic of the system (stability, settling time, noise 
sensitivity). 

According to Fig. 1, the frequency offset for the mth observation vector can be 
expressed as 

s 

z _j _. 2n&£\m-l\(i) 
= 2^ ln 732 g 7 2Mg H' 9) iJW^W* N ( 25 > 

where As[m-l]i(s[m-l]-s)[m-2])G and x m \i]±x[m{N +L)+S + i] The term inside the sum is 
a stochastic quantity and does not have a closed form expression. However, for 
reasonably high signal to noise ratio it can be well approximated by its expected value 

Therefore, the expression inside the sum can be written as 

^Ag[«-l]^ 



Z-l 2Ae[m-\\i) -J \— '— .2aA* m-lVi-l) Sin 

Z- 7, 1 — e " -j — - — - 

e N = — - , — ,- = e 2N 

,=0 l_ e J n sin 



V 2N j 



(Ae[m-l] 
2N 



27cM[m-\\L-l) L 



A£[m-l]->0 e 2N (26) 
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Substituting (26) into (25), after some algebraic manipulations, the frequency offset of 
m'th observation vector becomes 

e[m]=-(e[m-l]-s[m-2])^(S + ^) (27) 

The above equation represents a second order finite difference system in which its 
5 dynamic can be obtained from solving the following equation 
s[m]+fis[m - \]-pe[m - 2] = 0 

where fi- ^3 + ~-j . Clearly, the solution to the above finite difference equation is the 
form of 

e[m]= Cl {z,) m +c 2 {z 2 f (28) 

10 

where z l 2 = 2 are two dynamical modes of the system. The smaller root 

(negative) results in a high frequency oscillation in the frequency offset estimate. 
However, as is shown in a computer simulation, infra, this term is filtered out by the 
moving average filter. To assure stability, the gain block should be set such that both 
1 5 poles lie inside the unit circle. 



P max + V/? 2 max+4/?max|< 2 (29) 



where 



A Gf a L-V\ G(„ r L-l} 



(30) 



& N\ 2 ) N\ 2 
It is known that when the timing offset parameter is not known to the receiver or 
20 if the noise PDF differs from Gaussian distribution, finding the optimum estimator (ML, 
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MVU, CRLB) may not be an easy task. However there exists a moment estimator which 
provides a consistent estimate for estimation of frequency offset regardless of noise 
distribution and timing offset values. Although there is no optimum criterion associated 
with the moment estimator due to its simple structure, it is frequently used as an initial 
estimate for other estimators such as ML estimator. Consider a sequence of first N+L 
samples of vector x. Using the autocorrelation of Mi entry of this vector satisfies the 
following identity Moment Estimator 



a s e 



k£Q ( 31 ) 



Using Base Rule, the expected value of the above function (with respect to parameter k) 
can be expressed as 



r 

xx 



(32) 

= a } exp " J l7t€ pr (JceQyOpr (k$Q)= j~-ct } exp " J 2 ™ 



Substituting the Mh autocorrelation lag with its natural estimator, the moment estimator 
for frequency offset under uncertain timing offset can be found as 

s m om^\^{4\ (33) 
where the statistic Tg(x) is defined as 

t.w=^^V*M*+^1 (34) 

3 La} k=0 1 J L J 

Statistical assessment of moment estimator is a formidable task over entire range of SNR. 
However, for relatively high SNR, the random observation vector is heavily concentrated 
about its mean. Using the statistical linearization, and a first-order Taylor expression of 
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the estimator about its mean to obtain the variance of estimate. In doing so, one 
substitutes for the random variable x in (34) and obtains 



t 3 ([*])=/(M) = £{z*r;s*- W N 



k=0 



4k]) 



27t(k+N)s 
s*[k-S+N]exp J N 



J — ^ t r ' +w*\k+N~\ 



}(35) 



where the signal (s) and noise (w) vector are defined as 

s A[s[o]]-"j[jV + £-l]] 
w A[w[0]]...w[A^ + L-l]] 



(36,37) 



By virtue of the above equation, the expected value of observation vector for a 
fixed realization of signal vector s would be E[e mom ]= f(s) . Then performing a first 
order Taylor expansion of f([s.w]) about the point £w[ x ] yields 

smom =f([sM+Vwf([sM*w=f(s,0]) + L+ Z~ l ^&%w= 0)^ (38) 

n=Q M"} 1 

Taking the derivative of (35) with respect to wfnj and setting w=0, results in 



dh 



1 JV+X-1 n „ , 27te(n+N) 

w=0=^( I s[i-&]s*[i-& + N]e-J 27re T l s*[n + N-3]e * 

27U i-Q 



The second term in (38) represents the contribution of noise in the estimate. 

Knowing that noise samples are iid with power of a% , the variance of estimate can be 
obtained from 



varl 



(eh 



L+N-l 

z 

»— 0 



itJV-1 



r N+L-l 



(39) 
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For sufficiently large block lengths (AO, the above term can be well approximated as 



The resemblance between the estimators of equations (21) and (33) reveals a 
unified structure which characterizes the ML, MVU and moment estimators which can be 



of the proposed estimators. Moreover, it allows an investigation of the effect of symbol 
timing error in the estimation of carrier offset for each individual estimator. Comparing 
the MVU estimator given in (21) to the moment estimator in (33) reveals some 
similarities in the structure of the estimators. Clearly, both moment and MVU estimators 
10 use the same mapping function, namely the log function, to project the data statistics into 
the estimation domain. The only difference is in the form of statistics used for each 
scheme. 

FIG. 2 is a block diagram depiction of the ML estimator disclosed in the cross- 
referenced, copending application and the MVU and moment estimators of the instant 

15 invention. An autocorrelation operation on the observation vector x[k] is performed by 
blocks 30 and 34; as shown in FIG. 2, both estimators of the instant invention obtain the 
statistic by correlating the samples with the Mh delayed samples. This operation is 
performed by using a moving average (MA) filter (shown as a block 36) in the structure 
of estimators. However, the MVU and moment estimators use different upper and lower 

20 bound for the MA filter. In the moment estimator, the averaging is performed over the 
first N+L samples of the observation vector. This would remove the requirement of 
knowing the exact timing offset parameter in estimation of carrier frequency offset. 



var! 




(40) 



5 classified into a single unified structure. This provides a unique framework for analysis 
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However, the estimate obtained from using this estimator results in a less accurate 
estimate (more variance) in comparison to MVU estimate. On the other hand, the MVU 
estimator requires the knowledge of symbol timing in the estimation of carrier frequency 
offset, as shown in FIG. 2 by a block 38 conducting $ to the MA block 36. Finally, as in 
5 FIG. 1 , a block 40 and a block 42 perform the natural logarithm and imaginary, 

respectively, operations on the symbol stream before generating the estimate of carrier 
frequency offset s. 

Although the resemblance between MVU and ML estimators may not be as 
evident as that of MVU and moment estimator, it can be shown that ML estimator can 



10 also be classified into the same family. Knowing the fact that 3jlog7j^ ^ZT^Q) the 
ML estimator can be expressed as 



where the parameters a and J3 are functions of 3 ml and can be obtained from the 
expression for Ti(x,S), supra . Thus, the ML estimator falls into the same family of 

15 estimators. FIG. 2 also displays in tabular form these estimators, their respective 
summation upper and lower bounds, and their respective averaging intervals. It is 
noticed that the ML estimator provides the upper and lower bound of the moving average 
filter (36) by extracting the timing parameter from the likelihood function. Although the 
ML estimator has the advantage of exploiting the entire bandwidth by removing the 

20 requirement for having pilot tones, the symbol timing estimate obtained from the ML 

estimator has a larger confidence interval. This may result in a considerable performance 
degradation in comparison to the pilot-based schemes. 




(41) 
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Computer simulation is used to assess the performance of the three estimators for 
synchronization of an OFDM system. The variance of estimator is used as a performance 
measure through the study. The simulation parameters used are typical of the OFDM and 
digital audio broadcast (DAB) environments. More specifically, the chosen FFT size (N) 
5 for OFDM is 64. Unless otherwise specified, the length of cyclic prefix (L), signal to 
noise ratio, and frequency offset are set to 8, 20 dB and .01, respectively. Monte Carlo 
simulation is used to evaluate the performance of the three estimators. 

FIG. 3 graphically depicts a comparison between the performance of the ML 
frequency offset estimator disclosed in the related, copending application with the prior 

10 art ML frequency offset estimator given in the van de Beek, et al, reference over the 
range of timing offset parameter (& e [l, N + L] ) 

FIG. 4 graphically depicts the performance of the MVU frequency offset 
estimator of the instant invention under complete knowledge of timing offset error. A 
careful examination of the variances reveals that the gap between MVU estimator and 

15 CRLB tends to zero as SNR increases. Also as illustrated in FIG. 4, the departure from 
CRLB happens rapidly as SNR goes below a threshold. The threshold also depends on 
the length of CP and is moved toward lower SNRs as L increases. This can be justified in 
terms of having more observation samples in estimating the unknown parameter. The 
choice of cyclic prefix length L represents a tradeoff between data rate reduction and 

20 performance (lower variance). Increasing L brings the performance of MVU estimator 
closer to the CRLB, nevertheless, it could result in a considerable data rate reduction due 
to the redundancy introduced by CP. 
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FIG. 5 graphically depicts the closed loop performance of the MVU estimator of 
the instant invention by comparing the frequency offset estimate of the closed loop MVU 
estimator and the analytical derivation given in equation (28). It is clear that the 
simulation result very closely resembles the analytical model, thus consolidates the 
5 approximate model of the closed loop system. 

Numerous variations and modifications will become evident to those skilled in the 
art once the disclosure is fully appreciated. It is intended that the following claims be 
interpreted to embrace all such variations and modifications. 

10 WHAT IS CLAIMED IS: 

1. A method of estimating carrier frequency offset error in a received sample bit 
stream including an observation vector (OV), having an observed carrier frequency 
timing offset s, and a plurality of data-symbol frames, having a symbol timing offset 
error 3, comprising the steps of: 

15 generating a probability density function (PDF) based on said OV; and 

generating from said PDF an estimate of carrier frequency offset error, 
£mvu I being a minimum variance unbiased (MVU) estimator. 

2. The MVU carrier frequency offset error estimation method of claim 1 wherein 
20 said OV comprises an L-bit cyclic extension portion and a first and a second N-bit 

synchronization frame, and wherein said PDF comprises a first term, pi, based on said 
timing offset Q being within the span 1 to N and a second term, p2, based on said timing 
offset 0 being within the span N+l to N+L. 
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3. The MVU carrier frequency offset error estimation method of claim 2 wherein 
said received bit stream has uncorrected independent identically distributed random 
signal and noise sequence variables with power of a] and a\ , respectively, wherein 
5 said OV is denoted x, and wherein said MVU estimator smvu I a is the conditional 
expectation of a second moment estimator, said second moment estimator given by 




s k = 3 



4. The MVU carrier frequency offset error estimation method of claim 3 wherein 



1 0 said MVU estimator smvu I s is given by 



£ mvu 1,9 - E s 



Tfafi) J In 





where 3 is the imaginary operator and where 




8-N-l N+L-l 
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5. A method of synchronizing a received sample bit stream, comprising the steps 



of: 



transmitting at a transmitter said bit stream including an observation 



vector (OV); 
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receiving and sampling at a receiver said bit stream, said sampled bit 
stream including, said OV with an observed carrier frequency offset s, and a 
plurality of data-symbol frames, having a symbol timing offset error $ ; 

generating a probability density function (PDF) based on said OV; 

generating from said PDF an estimate of carrier frequency offset error, 
SmvuI being a minimum variance unbiased (MVU) estimator; and 

synchronizing said received bit stream by said MVU estimate of carrier 
frequency offset. 

6. The synchronization method of claim 5 wherein said OV comprises an L-bit 
cyclic extension portion and a first and a second N-bit synchronization frame, and 
wherein said PDF comprises a first term, pi , based on said observed timing offset S 
being within the span 1 to N and a second term, p2, based on said observed timing offset 
& being within the span N+l to N+L. 

7. The synchronization method of claim 6 wherein said received bit stream has 
uncorrelated independent identically distributed random signal and noise sequence 

variables with power of a] and a 2 n , respectively, wherein said OV is denoted x, and 

wherein said MVU estimator SmvuIo is the conditional expectation of a second moment 
estimator, said second moment estimator given by 
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8. The synchronization method of claim 7 wherein said MVU estimator smvu I a is 
given by 



£ mvu 1$ = E 



f 




€ 

I 






i 




In 



7i(x s v9) I 



5 where 3 is the imaginary operator and 



where Tj(x,.9)= 



.S4V-1 



A=0 



N+IA 

S 



9. A method of estimating carrier frequency offset error in a received sample bit 
stream including an observation vector (OV) 5 having an observed carrier frequency 

10 timing offset s, and a plurality of data-symbol frames, said OV comprising an L-bit cyclic 
extension portion and a first and a second N-bit synchronization frame, comprising the 
steps of: 

generating the expected value of the autocorrelation of the kth entry of 
said OV; and 

15 generating from said expected value an estimate of carrier frequency offset 

error, s m om 5 being a moment estimator. 

10. The moment carrier frequency offset error estimation method of claim 9 
wherein said OV is denoted x, and wherein said moment estimator e mom is given by 
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s mom- w~ 




where the statistic To(x) is defined as 




1 1. A method of synchronizing a received sample bit stream, comprising the steps 



transmitting at a transmitter said bit stream including an observation 
vector (OV); 

receiving and sampling at a receiver said bit stream, said sampled bit 
stream including, said OV with an observed carrier frequency offset s, and a 
plurality of data-symbol frames, having a symbol timing offset error 3 ; 

generating the expected value of the autocorrelation of the kth entry of 
said OV; 

generating from said expected value an estimate of carrier frequency offset 
error, s mom , being a moment estimator; and 

synchronizing said received bit stream by said moment estimate of carrier 
frequency offset. 

12. The synchronization method of claim 1 1 wherein said OV is denoted x, and 
wherein said moment estimator s mom is given by 



of: 
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where the statistic Tg(x) is defined as 

3 La} k=0 1 J L J - 

13. A method of estimating carrier frequency offset error in a received sample bit 
5 stream including an observation vector (OV), having an observed carrier frequency 

timing offset s, and a plurality of data-symbol frames, said OV comprising an L-bit cyclic 
extension portion and a first and a second N-bit synchronization frame, comprising the 
steps of: 

selecting one estimator method from a group of three: a maximum 
10 likelihood (ML) estimator, a minimum variance unbiased (MVU) estimator, a 

moment estimator, based on said OV; and 

generating from said selected estimator an estimate of carrier frequency 
offset error. 



15 14. The carrier frequency offset error estimation method of claim 13 wherein said 

MVU estimator method comprises the steps: 

generating a probability density function (PDF) based on said OV; and 
generating from said PDF an estimate of carrier frequency offset error, 
smvuU> being said MVU estimator. 

20 

15. The carrier frequency offset error estimation method of claim 14 wherein said 
OV comprises an L-bit cyclic extension portion and a first and a second N-bit 
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synchronization frame, and wherein said PDF comprises a first term, pi, based on said 
timing offset $ being within the span 1 to N and a second term, p2, based on said timing 
offset 9 being within the span N+l to N+L. 



10 



16. The carrier frequency offset error estimation method of claim 15 wherein said 
received bit stream has uncorrelated independent identically distributed random signal 
and noise sequence variables with power of a 2 s and a\ , respectively, wherein said OV is 
denoted x 9 and wherein said MVU estimator SmvuIs is the conditional expectation of a 
second moment estimator, said second moment estimator given by 



£ x[k]x*[k + N] 



2n {La} k = & 
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17. The carrier frequency offset error estimation method of claim 16 wherein said 
MVU estimator smvu I s is given by 





f 




MVU^ = E 


£ 

V 





3 

~2k 



r , s+l-\ 



2n 



Lai 



where 3 is the imaginary operator and where 

k=& 

s~m N+IA 



k=0 



k=3 
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18. The carrier frequency offset estimation method of claim 13 wherein said 
moment estimator method comprises the steps: 

generating the expected value of the autocorrelation of the kth entry of 
said OV; and 

5 generating from said expected value an estimate of carrier frequency offset 

error, £ mom? being said moment estimator. 

19. The carrier frequency offset error estimation method of claim 18 wherein said 
OV is denoted x, and wherein said moment estimator s mom is given by 



10 



8 mom- 




where the statistic To (x) is defined as 




15 



29 



LTI014 



ABSTRACT OF THE DISCLOSURE 
A class of non data-aided cyclic based robust estimators for frequency offset estimation 
of multi-carrier systems is disclosed. The use of sufficient statistics provides a minimum 
variance unbiased (MVU) estimate of the frequency offset under complete knowledge of timing 
offset error. The Neyman-Fisher factorization theorem and Rao-Blackwell-Lehmann-Scheffe 
theorem are used to identify the sufficient statistic and appropriate mapping functions. It is 
shown that there is but one function of the sufficient statistics which results in the minimum 
variance estimate among the possible class of cyclic-based estimators. Also, a moment estimator 
of frequency offset is provided to obtain a consistent estimate of carrier offset under uncertain 
symbol timing error. The moment estimator does not rely on any probabilistic assumptions. 
Thus, its performance is insensitive to the distribution of the additive noise. A unified structure 
characterizing both the MVU and moment estimators, as well as a maximum likelihood estimator 
of a related, copending application is disclosed. 
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